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We report the result from a search for bursts of gravitational waves using data 
collected by the cryogenic resonant detectors EXPLORER and NAUTILUS 
q \ during the year 2001, for a total measuring time of 90 days. With these data 

we repeated the coincidence search performed on the 1998 data (which showed 
a small coincidence excess) applying data analysis algorithms based on known 
', physical characteristics of the detectors. With the 2001 data a new interesting 

coincidence excess is found when the detectors are favorably oriented with 
respect to the Galactic Disk. 
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1. Introduction 



Cryogenic gravitational wave (GW) antennas entered into long term data taking operation in 1990 (EX- 
PLORER [p), in 1991 (ALLEGRO @|), in 1993 (NIOBE [|]), in 1994 (NAUTILUS [§]) and in 1997 
(AURIGA [||]), with gradual performance improvements over the years. 

Analysis of the data taken in coincidence among all cryogenic resonant detectors in operation 
during the years 1997 and 1998 was performed No coincidence excess was found above background 
using the event lists produced under the protocol of the International Gravitational Event Collaboration 
(IGEC), among the groups ALLEGRO, AURIGA, EXPLORER / NAUTILUS and NIOBE. 

Later [0], a coincidence search between the data of EXPLORER and NAUTILUS was carried out 
by introducing in the data analysis considerations based on physical characteristics of the detectors: the 
event energy and the directionality. The result was a small coincidence excess when the detectors were 
favorably oriented with respect to the Galactic Centre. 

Here we extend our analysis to new data obtained in the year 2001, when both EXPLORER and 
NAUTILUS were operating at their best sensitivity, using the same procedures applied for the previous 
analysis [f7|]. As previously done in ref. [^] we shall sometimes use the word probability, although 
we are well aware that its significance might be jeopardized by any possible data selection. With this 
proviso we shall use probability estimations in comparing different experimental conditions. 



Table 1: Main characteristics of the two detectors in the year 2001. The axes of the two detectors are aligned to within a few 
degrees of one other, the chance of coincidence detection thus being maximized. The pulse sensitivity for both detectors is of 
the order of h ~ 4 10~ 19 for 1 ms bursts. 
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2. Experimental data 

The resonant mass GW detectors NAUTILUS, operating at the INFN Frascati Laboratory, and EX- 
PLORER, operating at CERN, both consist of an Aluminium 2270 kg bar cooled to very low tempera- 
tures. A resonant transducer converts the mechanical oscillations into an electrical signal and is followed 
by a dcSQUID electronic amplifier. The bar and the resonant transducer form a coupled oscillator system 
with two resonant modes. 

With respect to the year 1998 the following changes were made in the set up of the detectors: 
the NAUTILUS detector operated at a thermodynamic temperature of 1.5 K instead of 0.14 K; the EX- 
PLORER detector was equipped with a new transducer providing a larger bandwidth and consequently 
enhanced sensitivity. The characteristics of the two detectors are given in the Table [[} 

The data, sampled at intervals of 12.8 ms for NAUTILUS and of 6.4 ms for EXPLORER, are fil- 
tered with an adaptive filter matched to delta-like signals for the detection of short bursts [[§]. This search 
for bursts is suitable for any transient GW which shows a nearly flat Fourier spectrum at the two resonant 
frequencies of each detector. The metric perturbation h can either be a millisecond pulse, a signal made 
by a few millisecond cycles, or a signal sweeping in frequency through the detector resonances. This 
search is therefore sensitive to different kinds of GW sources, such as a stellar gravitational collapse, the 
last stable orbits of an inspiraling neutron star or black hole binary, its merging and its final ringdown. 

Let x(t) be the filtered output of the detector. This quantity is normalized, using the detector 
calibration, such that its square gives the energy innovation of the oscillation for each sample, expressed 
in kelvin units. 

For well behaved noise due only to the thermal motion of the oscillators and to the electronic 
noise of the amplifier, the distribution of x(t) is normal with zero mean. Its variance (average value of 
the square of x(t)) is called effective temperature and is indicated with T e ff. The distribution of x(t) 
is 

f(x) = 1 e'^TJ (1) 

V 27rl eff 

In order to extract from the filtered data sequence events to be analyzed we set a threshold in terms of a 
critical ratio defined by 

CR = (2) 

where cr(|x|) is the standard deviation of \x\ and \x\ the moving average, computed over the preceeding 
ten minutes. 

The threshold is set at CR=6 in order to obtain, in the presence of thermal and electronic noise 
alone, a reasonable number of events per day (see ref.||7|]). This threshold corresponds to energy E t = 
19.5 T e ff. When |x| goes above the threshold, its time behaviour is considered until it falls back below 
the threshold for longer than three seconds. The maximum amplitude and its occurrence time define the 

event. 




The searched events are the ones that are due to a combination of a GW signal of energy E s and 
of the noise. The theoretical probability to detect a signal with a given signal to noise ratio R s = tM 3 -, 
in the presence of a well behaved Gaussian noise is [Q] 

probability (R s ) = / e 5 cosh{\f R e ■ R s )dR e (3) 

where i? e is the signal to noise ratio for the event and Rt = trr*- = 19.5 for the EXPLORER and 
NAUTILUS detectors. 

The behaviour of the integrand is shown in fig. |I[ This figure shows the spread of the event energy 
due to noise for a given R s of the signal. It shows that signals with R s = 20 (near the threshold) have a 
probability of about 50% not to be detected, and signals with R s = 30, rather larger than the threshold, 
still have a probability of near 15% not to be detected. The distinction between the two concepts, signal 
and event, is essential for our analysis, as discussed in ref [|7|]. 

Computation of the GW amplitude h from the energy signal E s requires a model for the sig- 
nal shape. A conventionally chosen shape is a short pulse lasting a time of r g , resulting (for optimal 
orientation, see later) in the relationship 



1 1 kE s 

where / is the resonance frequency, L and M the length and the mass of the bar and r g is conventionally 
assumed equal to 1 ms (for instance, for E s = 1 mK we have h = 2.5 10~ 19 ). 

3. Data selection 



All events which are in coincidence within a time window of ±5 s with signals observed by a seismome- 
ter are eliminated. This criterion cuts about 8% of the events. 




Fig. 2: The distributions of the hourly averages of T e // in kelvin units for EXPLORER and NAUTILUS. We accept only the 
time periods with hourly averages T e // < 10 mK. 



It is observed that the experimental data are affected by non gaussian noise which, in some cases, 
cannot be observed with any other auxiliary detector. Thus a strategy is needed to select periods during 
which the detectors operate in a satisfactory way, as discussed in the following paragraphs. 

To this end we consider the quantity T e ff, which we used in two ways. The first was to compute 
T e ff by averaging x 2 over one hour of continuous measurements (T e ft), the second to consider the T e /f 
averaged during the ten minutes preceeding each event. For the hourly averages we show the distribution 
in fig.||. On observing this figure we decided to consider for the search for coincidences only the time 
periods with hourly averages smaller than 10 mK. The distributions for the T e ff averaged over the ten 
minutes preceeding each event are shown in fig.|| for EXPLORER and NAUTILUS. It will be noticed 
that the number of events is larger for EXPLORER than for NAUTILUS. This depends on the bandwidth 
A/ which is larger for EXPLORER(see Table [[J). 

On observing these distributions we decided, as conservative a priori data selection, to make a 
cut and accept only the events for which the corresponding T e ff was below 7 mK (and there was no 
seismometer veto). This meant regarding the bump at 10 mK, for NAUTILUS, as due to extra noise. We 
recall that in the previous search [^] with the noisier 1998 data the cuts on T e ff were made at between 
25 and 100 mK. 

From our previous experience we had learned that the detectors operate in a more stationary way 
when the noise temperature remains low for longer periods of time, because this indicates a smaller 
contribution of extra noise. To make a quantitative check on this we classified the data stretches in various 
categories, according to the length of the continuous periods having hourly averages T e fj < 10 mK, 
obtaining the figures shown in Table ||. From this table we clearly see that the longer is the time period 
of continuous operation with low noise the smaller is the number of events associated with a noise 
T e ff > 7 mK. 

Finally, in fig.^] we show the distribution of the event energies selected according to T e f / < 10 mK 
and to T e jf < 7 mK, for each event, belonging to periods with duration > 1 hour. We notice that, in 
spite of our selection criteria, we still have several events with large energy, which indicates the presence 
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Fig. 3: The distributions of the ten minute averages of T e //, before each event, in kelvin units for the EXPLORER and 



NAUTILUS events (with hourly T" e // < 10 mK). We accept only the events with the ten minute average T e ff < 7 mK. 
Note that the EXPLORER events are more numerous than the NAUTILUS events because of the different bandwidth. 



Table 2: In the first column we indicate the minimum time length of continuous operation, in the second column the total com- 
mon time of measurement. The last columns indicate the number of events, only those during the common time of operation, 
with hourly T e f / < 10 mK and the percentage of the events having the ten minute average T e f / > 7 mK. 
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Fig. 4: The distributions of the EXPLORER and NAUTILUS event energies in kelvin units. 

of extra noise, in addition to the thermal and electronic ones. The only way to eliminate this noise is by 
means of the coincidence technique. 



4. Searching for coincidences 

For the search for coincidences it is important to establish the time window. Using simulated signals and 



real noise, we characterized [ ]10| ] the dispersion of the time of the event around the time when the signal 
is applied. 

The standard deviation of the time dispersion for a given detector is 



ad = const 



1 1 



(5) 



for delta signals, where A/ is the detector bandwidth and the const = 0.28 is determined, for the 
EXPLORER and NAUTILUS detectors, by means of simulation [[Tl|]. For a coincidence analysis with 
two detectors we have 

(6) 



o „ 



We decided to take as coincidence window w = ±3a w , which is the most recent choice of the IGEC 
collaboration. Since each event has its own a w the value of w will be different for each coincidence. 
Note that the value of w is almost entirely due to the NAUTILUS detector, since EXPLORER has a 
much larger bandwidth; it turns out to be of the order of \w\ ~ 0.5 s, about one half of that used in 
previous searches for coincidences. With the use of 3a w we also take into account the uncertainty of the 
detector bandwidth and of the simulation procedure. 

Analysis in a coincidence search consists of comparing the detected number of coincidences at 
zero time delay (±w) with the background, that is with coincidences occurring by chance. In order 
to measure the background due to the accidental coincidences, using a procedure adopted since the 



beginning of the gravitational wave experiments [12], one shift the time of occurrence of the events of 
one of the two detectors a number of times. We shifted 100 times in steps of At = 2 s (uncorrected 
data), from -100 s to +100 s. For each time shift we get a number of (shifted) coincidences. If the time 



shift is zero we get the number n c of observed coincidences. The accidental background is calculated 
from the average number of the n s hi t t shifted coincidences obtained from the one hundred time shifts 



n 



Y^lOO 

2^=1 



n 



ihift(j) 



100 



(7) 



This experimental procedure for evaluation of the background has the benefit of handling the problems 



arising when the distribution of the events is not stationary (see reference []13[]), although this is not the 
case with the present 2001 data. 



5. Energy filter 

It is clear that if a coincidence between the two detectors is due to the arrival of a GW burst we expect 
the energies of the two coincident events to be correlated, and we can disregard all coincidences whose 
corresponding event energies are very different, according to the considerations illustrated in fig.[]J Thus 
we can apply an energy filter with the aim of reducing the background. 

The procedure for application of such an energy filter was set in our previous search for coinci- 
dences [J7|]. We considered signals of various energies E s . For each coincidence found we calculated 
the R s for each of the above signal energy E s using the known values of the (local ) T e fj of the two 
events. We then verified whether the two R e , for the two events of that coincidence, fell within the inter- 
val R s ± AR S , such that the two limits R s — AR S and R s + AR S delimitate (see fig.[I|) an area of 68% 
(about one standard deviation for well behaved noise) for a given value of E s ; that is, we verified the 
compatibility of the two events. We followed the same procedure for the shifted coincidences in order to 
estimate the background after application of the energy filter. In this way we reduced, for the 1998 data, 
the average number of accidental coincidences from n = 223 ton = 51. 

This procedure is useful, in particular, if the two detectors have different sensitivity, as in the case 
of the 1998 data, and, consequently, the event thresholds and the event energies are also different. In the 
case of two detectors with comparable sensitivity, however, one could also consider to compare directly 
the energies of the coincident events. For the 2001 data, although in this year the two detectors had 
comparable sensitivities (see fig.||), we decided not to change the procedure used for the 1998 data. We 
considered GW signals of energy E s ,in a range covering the energies of our events, i.e. E s from 5 mK to 
1 K in steps of 5 mK, and accepted the coincidence (at zero delay or at a shifted time) if the two events 
fell within the above interval R s ± AR S . 



6. Sidereal time distribution 

In our previous search for coincidences [[7|] we took into consideration the non-isotropic response of the 
detector to a GW burst. We had reasoned that, since extragalactic GW signals should not be detected 
with the present detectors, possible sources should be located in our Galaxy, or in the Local Group. If 
any of these sources exist we should expect a more favorable condition of detection when the detectors 
are oriented with their axes perpendicular to the direction of the potential source, since the bar cross- 
section is proportional to sin 4 (8), where 9 is the angle between the detector axis and the direction of the 
line joining it with the source. We did find a small coincidence excess when angle 6 with respect to the 
Galactic Centre was larger than a certain value (see fig.3 of Ref.[|7|]). The inconvenience of this method 
is that, due to the poor statistics, the result has to be presented in an integral type graph, which makes it 
difficult to appreciate the real statistical significance of the data. Furthermore, hypotheses must be made 
on the location of the GW source. 

In the present search for coincidences we extend the previous analysis as follows. We still make 
use of the directional property of the antenna cross-section. As the Earth rotates around its axis, during 
the day the detector happens to be variably oriented with respect to a given source at an unknown location. 
Thus we expect the signal to be modulated during the day; more precisely the modulation is expected 
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Fig. 5: Result with events in the long time periods (> 12 hour) of continuous operation. The upper graph on the left shows 
the number of coincidences n c indicated with the * and the average number n of accidentals versus the sidereal hour. The lower 
graph on the left shows the Poisson probability to obtain a number of coincidences greater than or equal to n c . The two graphs 
on the right show the result using the solar time in hours. We remark that the data points refer to independent sets of events. 



to have a period of one sidereal day (with one or more maxima) (see references [[14], 15]), since the GW 
sources, if any, are certainly located far outside our Solar system. 

The principal, key analysis is carried out with the events in the time periods of at least twelve 
hours of continuous data taking (see Table |2|) to which the energy filter is applied. 

Twenty-four categories of events are considered, one per each sidereal hour, the sidereal time 
referred to a position and orientation halfway between EXPLORER and NAUTILUS (this determines 
the zero local sidereal time which is not essential for the following considerations). Each category 
includes coincidences totally independent from those in the other categories. For each category in fig. 
H we report the number n c of observed coincidences, the average number n of accidental coincidences 
obtained by using the time shifting procedure and, given n, the probability p that a number > n c of 
coincidences could have occurred by chance. For comparison we also show a histogram produced with 
the same procedure using solar hours. 

One notice a coincidence excess from sidereal hour 3 to sidereal hour 5, which appears to have 
some statistical significance, as the two largest excesses occur in two neighboring hours (the events in 
each hour are totally independent from those in a different hour). We have n c = 7 coincidences in this 
two-hour interval and n = 1.7. On the contrary, no significant coincidence excess appears at any solar 
hour. 

The accidental coincidences always have a Poissonian distribution. To check this, we have con- 
sidered for all the above events the accidental coincidences obtained with ten thousand trials, by time 
shifting from -10000 s to +10000 s in steps of two seconds. The distribution of the number of acciden- 
tal coincidences is shown in fig. [f| The agreement between experimental and expected distributions is 
excellent. 

We repeat the analysis for the events belonging to the larger set of continuous data taking lasting 
one hour or more (first line of Table^|). We obtain the result shown in fig.f^. We notice that in the sidereal 
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Fig. 6: The distribution of the number of accidental coincidences obtained with ten thousand trials, with average number 
n = 25.32. The stars indicate the experimental distribution, the continuous line the expected Poissonian distribution. 
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Fig. 7: Result with events in time periods > 1 hour of continuous operation. As in fig.| 
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Fig. 8: Correlation between the event energies of NAUTILUS with those of EXPLORER for the eight coincidences occurred 
in the sidereal hour interval 3 to 5, in time periods > 1 hour. The correlation coefficient is 0.96. No energy filter was applied. 



two-hour interval the number of coincident events has increased from seven to eight and the background 
has become n = 2.6. 

We proceed now in performing a test, to check if this result is compatible with simultaneous 
physical excitation of the two detectors of possible non-terrestrial origin. We compare the energies of 
the coincident events: if the events in EXPLORER and NAUTILUS are due to the same cause we expect 
their energies to be correlated. This test must be done without applying the energy filter. Using the events 
in the time periods with duration > 12 hours and without applying the energy filter we still get seven 
coincidences. If we consider the time periods with duration > 1 hour we get eight coincidences. The 
event energies are very strongly correlated, as shown in fig.^. We also studied the energy correlation of 
the events of the accidental coincidences and found no correlation. 

We also performed a coincidence data analysis when no energy filter at all was applied, obviously 
expecting a larger accidental background. The result is shown in fig.^. We find that the coincidence 
excess in the time interval 3 to 5 sidereal hours still shows up, although less clearly, as expected. 

A different way to present these data is shown in fig.[l(J This figure shows that at sidereal hours 
outside the interval 3 to 5 hours the event energies are not correlated. Also it shows the particular 
behaviour in the 3 to 5 hour interval, when the energies of all coincident events are correlated. 

The eight events in the 3 to 5 hour period are listed in the Table ^[ We have verified that, using the 
cosmic ray detector of NAUTILUS, these events are not due to cosmic ray showers. 



7. Comparison with the 1998 data 

The analysis presented here differs slightly (eg. for the use of the sidereal time) from that applied 
previously ^J§. We therefore present the 1998 data also in terms of the sidereal time. We must consider 
that the 1998 data are noisier than the 2001 data. In particular the EXPLORER data have a noise T e ff 
ten times larger than that of the 2001 data, whilst the NAUTILUS noise was of the same order. 

For the coincidence search we change the window from w = ±ls (the IGEC choice at that time) 
used in the paper [|7|] to the present w = ±3o~ w used here. The result is given in fig. 11. An effect similar 
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Fig. 9: As in fig-Q with no application of the energy filter. The two graphs on the left refer to events in time periods with 
duration > 12 hours. The two graphs on the right refer to events in time periods with duration > 1 hour. 
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Fig. 10: Ratios of the energies (nautilus/explorer) for the events in coincidence belonging to periods with duration > 1 hour 
(no application of the energy filter) versus the sidereal hour. 



Table 3: List of the coincident events at sidereal hours between 3 and 5. No energy filter has been applied. St is the time 
difference between the two coincident events and E is the event energy. 
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Fig. 11: 1998 data with energy filter. The upper graph on the left shows the number of coincidences n c indicated with the * 
and the average number n of accidentals versus the sidereal hour. The lower graph on the left shows the Poisson probability of 
obtaining a number of coincidences greater or equal to n c . The two graphs on the right show the result using the solar time in 
hours. 
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Fig. 12: 2001 data with "a posteriori" choices (see text). As in figj7| for the two graphs on the left (periods > 1 hour) and 
as in ngjijfor the two graphs on the right (periods > 12 hour). Again we remark that the data points in the figure at different 
sidereal hours are independent one from the other. 

to the one found for the 2001 data is noticed, although weaker than that obtained with the less noisy 2001 
data. 



8. Robustness of the statistical analysis 

In any data analysis care must be taken to avoid choosing procedures which favour particular results. 
Thus we considered the possibility that our result be biased, although involuntarily, by any such choices. 
In the present analysis all choices were made "a priori" and already published in the scientific literature; 
in particular the IGEC choice for the coincidence window of ±3 o~ w , the energy filter (see ref. ^ and 
[0]) and the threshold for definition of an event (see ref.[|6|]). 

Nevertheless, we tested whether our present choices were indeed, to accident, most apt to produce 
the coincidence excess. This proved not to be the case. 

To find the most favorable parameters, that is the threshold and the energy filter parameter, we 
considered only the coincidences in the sidereal hour interval 3 to 5 and minimized the probability of a 
coincidence excess by chance. We found that the most favorable threshold for the definition of event is 
at Rt = 20.5, instead of 19.5. The most favorable parameter for the energy filter is 50%, instead of 68%. 

We thought it interesting to report the result for these most favorable choices, R t = 20.5 for the 



threshold and 50% for the energy filter. The result is shown in fig. 12. We notice an indication that the 
coincidence excess might extend to the sidereal hour interval 3 to 6, including two more coincidences in 
the period 5 to 6 sidereal hours (see also figJT0|). We want to remark that this hour was not included in 
the optimization process. 

As far as the coincidence window is concerned, we found (a posteriori) that the best choice for 
having a coincidence excess in the 3 to 5 sidereal hour interval is ±3. 5cr w , and any coincidence window 
from ±2.5a w to ±4(7^ gives comparable results. 
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Fig. 13: The position of the Galactic Disk and the loci of point sources perpendicular to the detector axis at 1 and 13 sidereal 
hour in the right ascension-declination plane. At sidereal hour 4.3 the locus tends to coincide with the line of the Galactic Disk. 
The large circle indicates the location of the Galactic Centre. 



9. Discussion and conclusions 

The IGEC search for coincidences [^] was performed without applying the event energy algorithms 
based on the event amplitude and on the directional properties of the detectors; it gave a no coincidence 
excess. With new data taken in the year 2001 and with improved sensitivity we repeated the coincidence 
search with the detectors EXPLORER and NAUTILUS (no other detector was in operation during the 
year 2001), applying data analysis algorithms based on known physical characteristics of the detectors, 
namely energy of the events and directionality of the detectors. We obtained a coincidence excess at 
sidereal hours between 3 and 5. 

At a given sidereal time the intersection of the celestial sphere with the plane perpendicular to the 



detector axis is a circle. We show in fig.|13| two of these circles (i.e., at 1 and at 13 sidereal hours) in 
the right ascension-declination plane. In the representation of figjl^ the line which indicates the point 
sources perpendicular to the detector axis (we call this the line of maximum sensitivity) moves to 
the right with sidereal time. This line intersects the location of the Galactic Centre twice a day (at 4.3 
and at 13.6 sidereal hours). Only once per day (at 4.3 sidereal hour) this line overlaps with the entire 
Galactic Disk. The overlapping takes place because of the particular orientation of the detectors on the 
Earth' surface (see Table [[]), and would not occur for a different azimuth angle of detector orientation. 

If all GW sources were concentrated in the Galactic Centre we would have found a coincidence 
excess twice a day. The coincidence excess occurs only once per day, just when the line of maximum 
sensitivity of the detectors overlaps with the Galactic Disk, as if GW sources were distributed in the 
Galactic Disk and not just located in its Center. 

As for the energy balance, in the year 2001 we find in the interval from 3 to 5 sidereal hour a 
coincidence excess of, very roughly, n c — n ~ 6 coincidences occurring in five days (two sidereal hours 
out of twentyfour, in a total time period of 1490 hours ~60 days). In terms of energy conversion into 
GW we have, very roughly, about one coincidence per day with a signal energy of about 100 mK. This 
corresponds, using the classical cross-section, to a conventional burst with amplitude h ~ 2 10 -18 and 



to the isotropic conversion into GW energy of 0.004 solar masses, with sources located at distance of 8 
kpc. The observed rate is much larger than the models today available predict, for galactic sources. We 
note, however, that our rate of events is within the upper limit determined by IGEC [0] for short GW 



bursts and by the 40m-LIGO prototype interferometer Q16Q for coalescing binary sources in the Galaxy. 

We think it is unlikely that the observed coincidence excess be due to noise fluctuations, but we 
prefer to take a conservative position and wait for a stronger confirmation of our result, before reaching 
any definite conclusion and claim that gravitational waves have been observed. Furthermore, although 
we have excluded that the events are due to cosmic ray showers (see section 6), we cannot completely 
rule out that they be due to some other exotic, still unknown, phenomenon. A possible way to distinguish 
GW from other causes is to measure other vibrational modes of the detectors and verify that, as predicted 
by General Relativity, only quadrupole modes are excited. This requires multimode detection (with bars 
or spheres) and to improve the signal-to-noise ratio of the apparatuses. 

We expect to collect new data with EXPLORER and NAUTILUS with improved sensitivity. We 
plan to repeat the same analysis with these new data and also with any other new data provided by other 
GW groups, those which operate the resonant detectors and those which operate or are beginning to 
operate the interferometric detectors GEO, LIGO, TAMA and VIRGO. 
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